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Processes involving flavor changing neutral currents (FCNC) provide excellent signatures with 
which to search for evidence of new physics. They have very small branching fractions in the 
Standard Model since they are highly suppressed by Glashow-Iliopoulos-Maiani (GIM) mecha- 
nism. They occur only through higher order diagrams, and new particles contributions can provide 
a significant enhancements, which would be an uniquevocal signs of physics beyond the Standard 
Model. In this paper we present the most recent measurements on FCNC processes performed by 
CDF and D0 Collaborations, while last section is devote to the charm physics at CDF. 
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1. Search for rare — > decay modes 

The FCNC decays IJ.^IJ.^ ^ occur in the Standard Model (SM) only through 

higher order diagrams and are further suppressed by the helicity factor, (m^/m^)^. The de- 
cay is also suppressed with respect to the B*' decay by the ratio of CKM elements, iVid/V^sl^. 
The SM expectations for these branching fractions are ^(6^ — )• /i^/^") = (3.42 it 0.54) x 10"^ 
and ^{B*^ — )• jU+jU") = (1.00ib0.14) x 10"^° [Q], which are one order of magnitude smaller than 
current experimental sensitivity. Enhancements to jX^pL^ occur in many new physics 

models [^. In the absence of an observation, limits on ^(B^ — )■ fX^fX^) are complementary to 
those provided by other experimental measurements, and together would significantly constrain 
the allowed supersymmetric parameter space. In general, the search for these rare decays is central 
to exploring a large class of new physics models. In the following we will present recent results 
on the search of these very rare decays from CDF |5]] and D0 [^, ^] Collaborations. Since two 
analyses are very similar, we will shortly describe only CDF analysis. Details on CDF and D0 
detectors can be found in 

At CDF in the offline analysis, the trigger selection [0, |]] is refined by applying a series of 
"baseline" requirements that substantially reduce the backgrounds while preserving the majority 
of the signal. We select two oppositely charged muon candidates within a dimuon invariant mass 
window of 4.669 < m^^ < 5.969GeV/c^ around the and B^ masses. Backgrounds from hadrons 
misidentified as muons are suppressed by selecting muon candidates using a likelihood function. 
This function tests the consistency of electromagnetic and hadronic energy with that expected for 
a minimum ionizing particle and the differences between extrapolated track trajectories and muon 
system hits [^|. In addition, backgrounds from kaons that penetrate through the calorimeter to the 
muon system or decay in flight outside the drift chamber are further suppressed by a loose selec- 
tion based on the measurement of the ionization per unit path length, dE/dx [^, [T^. The inputs 
to the muon likelihood and the dE / dx performance are calibrated using samples of7/i/A— )• pL^pL^ , 
D" — )• K^n^ and A — )• pn^ decays. To reduce combinatorial backgrounds the muon candidates 
are required to have transverse momentum relative to the beam direction pr > 2.0 GeV/c, and 
\Py^\ > 4 GeV/c, where p^^ is the transverse component of the sum of the muon momentum 
vectors. The remaining pairs of muon tracks are fit under the constraint that they come from the 
same three-dimensional (3D) space point. To achieve further separation of signal from background, 
we employ additional discriminating variables. These include the measured proper decay time. A; 
the proper decay time divided by the estimated uncertainty, Xjox, the 3D opening angle between 
vectors p^^ and the displacement vector between the primary vertex and the dimuon vertex, A0; 



and the B-candidate track isolation, / 1 11]. We require that A/a;^ > 2, A0 < 0.7 rad, and / > 0.50. 
The baseline selection reduces combinatorial backgrounds by a factor of 300 while keeping ap- 
proximately 50% of the signal events that are within the acceptance (geometric and kinematic 
requirements) of the trigger. A sample of — )• J /\\rK^ events is collected to serve as a normal- 
ization mode using the same baseline requirements, but including a requirement oi pT > I GeV/c 
for the kaon candidate and constructing the 5+ — >• J/xj/K^ — >• l^^l^ vertex using only the muon 
candidate tracks. 

For the final event selection we use the following discriminating variables: m^^. A, ^/o^, 
A0, /, \pj^\, and the pj of the lower momentum muon candidate. To enhance signal and back- 
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Figure 1: Distributions of Vn for simulated l^^li-^ signal and observed sideband events (on the left). 
The jJ-^IJ.^ invariant mass distribution for events satisfying all selection criteria for the final three ranges of 
VV (on the right). 



ground separation we construct a NN discriminant, Vn, based on all the discriminating variables 
except m^^, which is used to define signal and sideband background regions. The NN is trained 
using background events sampled from the sideband regions and signal events generated with a 
simulation described below. The Vn distributions of signal and sideband background events are 
shown in fig. [T[ 

For measuring efficiencies, estimating backgrounds, and optimizing the analysis, samples of 
— > jX^jX^ , 5+ — ^ J/yK^, and B — ^ h~^h~ (where are or K^) are generated with the 
PYTHIA simulation program [ |l2| ] and a CDF II detector simulation. The B-hadron pj spectrum 
and the / distribution of the B-hadrons are weighted to match distributions measured in samples of 
B+ J/\i/K+ and B° decays. 

We use a relative normalization to determine the B*' — H^jJ.^ branching fraction: 

as Es Sn Js 

where Ng is the number of — )• H^jJ.^ candidate events. We observe about ^ 19,700 — 
J/yK+ candidates. We use ^{B+) = ^{B+ J/yK+ }X+IX-K+) = (5.94 ±0.21) x 10"^ [111 



and the ratio of B-hadron production fractions /„//« = 3.86 ±0.59 Jlsj]. The parameter as («+) 
is the acceptance of the trigger and Es (£+) is the efficiency of the reconstruction requirements 
for the signal (normalization) mode. The reconstruction efficiency includes trigger, track, muon, 
and basehne selection efficiencies. The NN efficiency, En, only apphes to the signal mode. The 
expression for SS{B^ — )• jx^ ji^) is derived by replacing B° with B^ and the fragmentation ratio with 

fu/fd = 1- 

The ratio of acceptances a+ / a, is measured using simulated events, and its uncertainty in- 
cludes contributions from systematic variations of the modeling of the B-hadron pT distributions, 
the longitudinal beam profile, and from the statistics of the simulated event samples. The ratio 
of reconstruction efficiencies £+ / is measured using simultaneously data and simulation. Muon 
reconstruction efficiencies are estimated as a function of muon pj using observed event samples 
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of inclusive J/y ^ P^^P^ decays. Systematic uncertainties in this ratio largely cancel with the 
exception of the kaon efficiency from the B+ decay. The uncertainty is dominated by kinematic 
differences between inclusive J/y ^ M^M ^i^d £['(6") — )• IJ^^IJ^^ decays. The efficiency, En, is 
estimated from the simulation. We assign a relative systematic uncertainty on of 6% based on 
comparisons of NN performance in simulated and observed — J/yK^ event samples and the 
statistical uncertainty on studies of the B° pj and / distributions from observed — J/y'P decays. 

The expected background is obtained by summing contributions from the combinatorial con- 
tinuum and from B — )• h^h^ decays, which peak in the and invariant mass signal region and 
do not occur in the sidebands. The contribution from other heavy-flavor decays is negligible. We 
estimate the combinatorial background by linearly extrapolating from the sideband region to the 
signal region. The B — )• contributions are about a factor of ten smaller than the combina- 

torial background and are estimated using efficiencies taken from the simulation, probabilities of 
misidentifying hadrons as muons measured in a — )- nK data sample, and normalizations derived 
from branching fractions from refs. |^, 13]. The two-body invariant mass distribution of the sim- 



ulated B — )• candidates is calculated from the momentum of the hadrons assuming the muon 
mass hypothesis. The background estimates are cross-checked using three independent control 
samples: jJ-'^jJ.'^ events, ix^fx^ events with A < 0, and a misidentified muon-enhanced fX^IX^ sam- 
ple in which we require one muon candidate to fail the muon quality requirements. We compare 
the predicted and observed number of events in these samples for a wide range of Va? requirements 
and observe no significant discrepancies. 

The pL^li^ invariant mass distributions for the three different ranges are shown in fig. |l]. 
The observed event rates are consistent with SM background expectations. Using a data sample 
of 3.7 fb^' of integrated luminosity, collected by CDF II experiment, we extract 95% (90%) C.L. 
limits of ^(S? ^ At+M") < 4.3 X 10-*^ (3.6 x 10"^) [|] wd ^{B^ ^ pL+ pL') <1 .6x\Q-'^ (6.0 x 
10-'') [g, which are currently the world's best upper hmits for both processes. Assuming A^v = 1, 



from eq. LI we compute the single event sensitivity (SES) for the B° — )• H^jJ.^ decay mode as 
3.2 X 10"^ for all mass and Vn bins. The SES for the B^ — > M^M" is smaller than the expected SM 
branching fraction and we expect 1.2 events from — pL^pi^ decays with 0.7 events occurring in 
the highest sensitivity Vn bin. The result for 5^ — pi^pL^ is slightly in excess of the expected limit, 
which is ^(S^ M+M") < 3.3 x 10^^(2.7 x 10"**) at 95(90)% C.L., and this was estimated with 
no signal hypothesis. We calculate the P-value of the excess as 23% corresponding to 0.73a. 

D0 performs a similar analysis. Latest results are based on a data sample of integrated lumi- 
nosity of 2 fb" ' . The observed event rates are, also in this case, consistent with SM background 
expectations. We extract 95(90)% C.L. limit of ^(B° pi^pi^) < 9.3 x 10"^ (7.5 x 10"^) ||]. 
With 5 fb"^ of data collected by the D0 experiment, we have studied the sensitivity to the branch- 
ing fraction of B° — )• fX^fX^ decays. An expected upper limit on the branching fraction is =^(5° — )• 
pi+pi-) < 5.3(4.3) X 10-^ at the 95(90)% C.L. ^. 

We observe no evidence for new physics and set limits that are the most stringent to date, 
improving the previous results [^, ^, |l^] by a factor of 1.5 or more. These limits place further 
constraints on new physics models, and complement its direct searches. We expect the analysis 
sensitivity to continue to improve as we include larger data sets. 
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Figure 2: The B invariant mass of B+ /i:+^+/i"(on the left), (center) and B° 

(j) /X+jU^(on the right) on a data sample of 4.4 fb^^ of integrated luminosity. Solid, dashed, and dotted line 
shows total fit, signal PDF and background PDF, respectively, while dots with error bars are data. 



2. b ^ s t^l transition at the Tevatron 

The b ^ s i^i^ transition is a FCNC process, which, in the SM, proceeds at lowest order 
via either a Z/y penguin diagram or a W+W^ box diagram, as the — )■ M^M^ decays. The 
effective Wilson coefficients Cy, C9, and Cio describe the amplitudes from the electromagnetic pen- 
guin, the vector electroweak, and the axial-vector electroweak contributions, respectively. These 



amplitudes may interfere with the contributions from non-SM particles [15], therefore the transi- 
tion can probe the presence of yet unobserved particles and processes. More specifically, the lepton 
forward-backward asymmetry (Afs) and the differential branching fraction as functions of dilep- 
ton invariant mass (Mu) in the decays B — )• K*£^i^ differ from the SM expectations in various 
extended models []T6[]. The former is largely insensitive to the theoretical uncertainties of the form 
factors describing the decay, and can hence provide a stringent experimental test of the SM. The 
latter has been already determined [|^, |l^], however the precision is still limited by statistics. It 
can be used to extract the information on the coefficients associated with the theoretical models as 
well. For more details see ref. [|T9|]. 

Although the branching ratio of each b ^ s jJ-^fJ.^ decay is quite small as (f?(10^^), the 
decay is experimentally clean due to opposite sign muons. Among many b ^ s H^jJ.^ decays, 
the exclusive channels B+ — )• K^pL^ji and B" — )• K*^ii^pL'' have been observed and studied at 



Belle rtnl] and BaBar [IT^. However, the analogous decays S° — )• ^ jx^fi. and A° — > A 



have not been observed despite searches by CDF p(3| ] and D0 [21]. Electrons can be reconstructed 
at the TeVatron, but they involve additional experimental difficulties and they will be treated in next 
iterations of the analysis. 

CDF is currently updating the analysis on a data sample of 4.4 fb^^ of integrated luminos- 



ity [||]. We obtain a signal yield for the S° K*°ix+n^ decay mode equal to 98 ± 12 (101 
expected), with a statistical significance of about 10a and a signal yield for the B+ — > K^n^jj.^ 
equal to 119 it 16 (140 expected), with a statistical significance of about 9a. Significances are 
deterimined from the likelihood ratio to a null signal hypothesis. Search of the B° — )• M^M^ 
is still blind as shown in fig. 0. We expect a significant signal yield, also in this case, of about 31 
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events. Measurement of ^(M|^) and A^g is underway [23]. 

3. Search for rare — > decay mode 

The FCNC decay D" — )• pL^pi^ is highly suppressed in the SM by the nearly exact Glashow- 



Iliopoulos-Maiani cancellation. SM expects the branching fraction to be about 10 [24] from 
short-distance processes, increasing to about 4 x 10^ [24] by including long-distance processes. 



The prediction is many orders of magnitude beyond the reach of the present generation of experi- 
ments. The best published upper bound is 1.3 x 10^^ at the 90% C.L. from BaBar [25], while the 



world's best upper bound is 1.4 x 10"'^ at the 90% C.L. from Belle [^. 

However, new physics contributions can significantly enhance the branching ratio [24]. Some 
of these new scenarios could enhance the branching fraction to the range of 10"^ to 10~^° , and in 
particular R-parity violating SUSY could lift the branching fraction up to the level of the existing 
experimental bound. Similar enhancements can occur in K and B decays, but charm decays provide 



a unique laboratory to search for new physics couplings in the up-quark sector. Ref. [ ]27| ] shows 



that, in some scenarios, new physics contributions to — d" mixing can dominate, but yield a 
result consistent with the SM, the same new physics can contribute to — )• pi^jx^ , yielding a 
visible signal. 

Using a data sample of 360 pb^^ of integrated luminosity CDF searches for D'' — )• 
decays. A displaced-track trigger selects long-lived D" candidates in the pL^pL^ , tt+tt^, and K^n^ 
decay modes. We use the kinematically similar D" — )• TT+TT^ channel for normalization, and the 
Cabibbo-favored D" — )■ K 7r+ channel to optimize the selection criteria in an unbiased manner. 
We set an upper limit on the branching fraction SS{D^ — )- pL^pL^) < 3.0 x 10^^(2.1 x 10^^) at 
the 95(90%) C.L. [^]. CDF is currently analyzing a data sample of about 5 fb^^ of data (already 
on tape), and it will integrate 8(10) fb^^ by the end of 2010(2011). We expect to significantly 
approach, by the end of Run II, the interesting region 10^^ for the D" — p^pi^ decay. 

4. Charm IMixing 

In the SM, the decay D" — K^K^ proceeds through a doubly Cabibbo-suppressed (DCS) 
"tree" diagram, and may also result from a mixing process (D° -H- D^), followed by a Cabibbo- 
favored (CF) decay (D*' — K^n^). The DCS decay rate depends on Cabibbo-Kobayashi-Maskawa 
quark-mixing matrix elements and on the magnitude of SU(3) flavor symmetry violation [^]. 
Mixing may occur through two distinct types of second-order weak processes. In the first, the 

evolves into a virtual ("long-range") intermediate state such as 71+71^, which subsequently 


evolves to a Z) . The magnitude of the amplitude for long-range mixing has been estimated using 



strong interaction models [|30|], but has not been determined using a QCD calculation from first 
principles. The second type of second-order weak process is short-range [|1]], with either a "box" 
or "penguin" topology. Short-range mixing is negligible in the standard model. However, exotic 
weakly interacting particles could enhance the short-range mixing and provide a signature of new 



physics [32] 



The ratio /? of D" — )• K^n to D" K ti^ decay rates can be approximated as a simple 
quadratic function oft /x, where t is the proper decay time and T is the mean D° lifetime. This form 
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is valid assuming CP conservation and small values for the parameters x = AM /F and y = AT /IF, 
where AM is the mass difference between the meson weak eigenstates, AF is the decay width 
difference, and F is the average decay width of the eigenstates. Under the assumptions stated above, 

R{t/T) = Ro + VR^y' {t/T) + {t/T)\ (4.1) 

where Rd is the squared modulus of the ratio of DCS to CF amplitudes. The parameters x' and y' 
are linear combinations of x and y according to the relations 

x' = X cos 5 +3^ sin 5 and / = — x sin 5 +3^ cos 5, 

where 8 is the strong interaction phase difference between the DCS and CF amplitudes. In the 
absence of mixing, x! = / = and R{t /x) = Ro- 

We use a signal of 12.7 x 10^ D" — K^n^ decays with proper decay times between 0.75 and 
10 mean lifetimes. The data sample was recorded with the CDF 11 detector and corresponds 
to an integrated luminosity of 1.5 fb^'. We search for — D mixing and measure the mixing 
parameters to he Rd = (3.04 ±0.55) x lO^^,/ = (8.5 ±7.6) x 10"^ and x'^ = (-0.12 ±0.35) x 



10 [33]. Bayesian probability contours in the x — / plane shows that the data are inconsis- 



tent with the no-mixing hypothesis with a probability equivalent to 3.8a 1330. This result is in 



agreement with the current B -factories measurements []34{], and it has a comparable sensitivity. 

Using the current data sample, corresponding to 4.4 fb^' of integrated luminosity, we recon- 
struct about 24 X 10^ D° — )• K^n^ decays. This is approximately a factor 2 better in statistics with 



respect to the published results [p3[]. Assuming no analysis improvements and assuming as central 
values those published in [33], we expects to reach 5a level of significance and to provide a precise 
measurement of mixing parameters x'^ and y' at the next iteration of this analysis. CDF is currently 
taking data with the goal of integrating 8 fb^' by the end of 2010. If the extension of Run 11 will 
be approved, CDF will collect on tape about 10 fb^^ by the end of 2011, of good data for physics 
analysis. This scenario is very promising considering that CDF has already the world's largest data 
sample, and it is taking new charm data with an unprecedent rate. 
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